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Abstract 

We use boson-fermion mapping to show that the single-particle momentum dis- 
tribution in a one-dimensional gas of hard point-like bosons (Tonks gas) inside a 
harmonic trap decays as p~ 4 at large momentum p. The relevant integrals expressing 
the one-body density matrix are evaluated for small numbers of particles in a simple 
Monte Carlo approach to test the extent of the asymptotic law and to illustrate the 
slow decay of correlations between the matter-wave field at different points. 
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1 Introduction 

The realization of Bose-Einstein condensates of lower dimensionality in suit- 
able optical and magnetic traps [1,2], and advances in atom waveguide tech- 
nology with potential applications to atom interferometry and integrated atom 
optics [3-7] , give special motivation to theoretical studies of dilute atomic flu- 
ids in a regime where the quantum dynamics becomes quasi-one-dimensional 
(ID). The condensate fraction is depleted in this situation even at zero temper- 
ature, and the standard stationary or time-dependent Gross-Pitaevskii theory 
loses its usefulness [8-10]. An extreme limiting example is provided by the 
so-called Tonks gas of hard point-like bosons, in which the transverse con- 
finement is taken to be so strong that the dynamics reduces to strictly ID 
motion. The occupation of the lowest orbital behaves in this case with the 
number N of particles as N n with n < 1 [11]. The momentum distribution 
still has a strong peak on approaching zero momentum [8,11], and coherence 
effects have been shown to be present in the form of Talbot recurrences follow- 
ing an optical lattice pulse [12] and of dark soliton-like behaviour in response 
to a phase-imprinting pulse [13]. 
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A number of properties of the Tonks gas of impenetrable bosons can be cal- 
culated exactly through a mapping into an ideal ID gas of spin-polarized 
fermions [14]. In both fluids each particle behaves as a hard point which is at 
any time free to move on a segment comprised by its two immediate neighbours 
and can exchange with either of them on contact. Considering in particular 
the many-body ground state, the boson-fermion mapping consists of the iden- 
tity between the bosonic wave function and the modulus of the fermionic wave 
function [14]. This of course implies that, whereas all configurational probabil- 
ity distributions are the same in the two fluids, their single-particle momentum 
distributions are in general completely different [15]. 

A visualization of the ground state can be obtained by varying the position 
of one particle while all others are kept fixed. The fermionic wave function 
varies smoothly from positive to negative across each node so as to mimimize 
the kinetic energy, and hence the bosonic wave function, being everywhere 
non-negative, has a cusp at each node. We show in this paper that from this 
cusp condition the momentum distribution of the Tonks gas under harmonic 
confinement acquires a long tail at high momentum p, decaying asymptotically 
as p~ 4 . We also gauge the extent of this asymptotic behaviour by Monte Carlo 
calculations of the momentum distribution for the case of two particles in a 
harmonic well, where high statistical accuracy and simple analytical results 
are easily obtained. 



2 Asymptotic tails of the bosonic momentum distribution 



The single-particle momentum distribution n(p) at momentum p = hk is the 
Fourier transform of the one-body density matrix p(x, x') with respect to the 
relative coordinate r = x — x', averaged over the centre-of-mass coordinate 
R = (x + x')/2. The ground-state wave function for a Tonks gas of N bosons 
under longitudinal harmonic confinement has been calculated by Girardeau et 
al. [11], with the result 

N 

^o(x 1 ,...,x N ) = C N n \x k -x j \Y[exp(-x 2 j /2x 2 ho ). (1) 

l<j<k<N i=l 



Here, Cn is a normalization constant and Xho = (h/muj) 1 / 2 is the harmonic 
oscillator length determined by the trap frequency u and the particle mass m. 
Taking x ho as the unit of length, the momentum distribution is then given by 

+oo +oo +oo +oo 

n(p) = M N J dR J dr exp[-ikr - (R 2 + r 2 /4)] J dQ 2 ■ ■ ■ J dQ N 
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where Mn is a factor depending on the number of particles. 

Let us first treat the case N = 2. The momentum distribution is given by 



-oo +oo 



nip) = M 2 J dR J dQ 2 exp{-R 2 - Q 2 2 )I{Q 2 , R; k) (3) 
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where 

+oo 



I{Q 2 ,R;k)= J dr\{Q 2 -R) 2 -r 2 /A\exp{-tkr-r 2 /A). (4) 



The integrand has two singularities of the type \r ± a\ and we can use a 
theorem given in the book of Lighthill [16] to evaluate its asymptotic value in 
the limit k — > oo. Using the result / dz\z \ exp(ikz) = —2/k 2 we get 

lim I{Q 2 , R; k) = _ 7 \9l_l3. cos [2A;(Q 2 - R)\ exp[-(Q 2 - R) 2 }. (5) 
fc^oo /e z 



Introducing the notations X = (Q 2 — R)/ \[2 and Y = (Q 2 + R)/ y/2, the 
asymptotic behaviour of the momentum distribution is then given by 

lim n{p) = Z -^-N 2 J dYexp(-Y 2 ) J dX\X\ cos(2 v / 2A;X) exp(-3X 2 )(6) 



and Lighthill's theorem can again be employed to evaluate the integral over 
X. The final result is 

W-m nip ) = 2\ -± 4 , 7) 



having used the value of the normalization constant C 2 reported by Girardeau 
et al. [11]. 

In applying the above procedure to the calculation of the asymptotic value of 
Eq. (2) for arbitrary N, we notice that the same pair of singularities arises 
for each value of the index i and that they all contribute equally to the final 
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result. Choosing i = 2, we thus find 

lim n(p) = —7- (8) 



where the coefficient A N is determined by a multiple integral, 



+00 +00 +00 / N \ 

A N oc — (A" — 1) A/jv I rfFexp(-F 2 ) / dQ 3 . . . J dQ N expl-^Q 2 j 



X 



II (Qk - Qjf (Qk - Y/V2) (Qj - Y/V2) 

3<j<k<N 



(9) 



We have thus shown that the momentum distribution of the Tonks gas in 
a harmonic trap decays asymptotically as the inverse fourth power of the 
momentum. 



3 Numerical results by a Monte Carlo method 



In this section we evaluate the one-body density matrix and the momentum 
distribution function for a small number of particles by a simple Monte Carlo 
method. The density matrix involves the integrals over (Q2, ■ ■ ■ ,Qn) enter- 
ing the second line in Eq. (2). These are calculated by (i) extracting for each 
variable Qi a random number distributed along a Gaussian of unitary width, 
and (ii) averaging over all trials the corresponding values of the factor con- 
tained in the curly brackets in the integrand. The two further integrals that 
are required for the evaluation of n{p) are calculated by a similar procedure. 
A more refined numerical method would be needed to obtain high statistical 
accuracy for a number of particles N > h. 

We report in the following our numerical results for N = 2. Figure 1 shows the 
function p(x, 0) as a function of x for a pair of hard bosons, after suppression 
of the Gaussian decay imposed by the confinement through multiplication 
by [p(x)p(0)]~ 1/ ' 2 where p(x) is the particle density. The comparison with the 
same function for a pair of non-interacting fermions, for which a direct analytic 
calculation yields p(x, 0)/[p(x)p(0)} 1 ^ 2 = (1 + 2x 2 )~ 1 ^ 2 , illustrates the role of 
statistics in determining the rate of decay of the density matrix. It may be 
remarked that the density matrix for the Bose gas reflects its spatial coherence 
and can be determined experimentally by measuring the interference pattern 
of two matter-wave beams emitted from two spatially separated regions of the 
trap [17]. 
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Figure 2 reports the corresponding results that we obtain for the momentum 
distribution n(p) in the case of two Bose or Fermi particles. The analytical 
result for fermions, with n(p) described by the function n~ 1 / 2 (l+2y 2 ) exp(—y 2 ) 
where y = p/ '(hmuo) 1 / 2 , is shown by the full line and serves to illustrate the 
accuracy of the Monte Carlo calculation. 

Finally, Figure 3 tests the range of applicability of the asymptotic law given 
in Eq. (7) by reporting it in a log-log plot against the momentum distribution 
obtained in two Monte Carlo runs for a pair of hard bosons. 



4 Concluding remarks 

In this work we have demonstrated two aspects of the correlations which arise 
from Bose statistics in a ID gas of hard, point-like bosons inside a harmonic 
trap. A very remarkable feature of this system is the extremely slow decay of 
the high-momentum tail in the single-particle momentum distribution, which 
we have determined analytically and verified numerically. We have also shown 
that the one-body density matrix decays very slowly, after elimination of the 
Gaussian factors that are imposed by the trap. 

A consequence of the p~ 4 tail in the momentum distribution is that the mean 
square fluctuation of the kinetic energy of the Tonks gas, which involves the 
fourth moment of n(p), diverges under harmonic confinement. This property 
is even more remarkable if one considers that the kinetic energy density of 
the Tonks gas is the same as that of its fermionic map. This latter result, 
which can be demonstrated analytically in full detail from the bosonic and 
fermionic ground-state wave functions [18], may be viewed most directly to be 
a consequence of the basic theorem of density functional theory: the kinetic 
energy density is a unique functional of the particle density, which remains 
unchanged in the boson-fermion mapping. 

We acknowledge support from INFM through PRA2001. 
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Fig. 1. Monte Carlo results for the one-body density matrix p(x, x')/[p(x)p(x')] 1 ^ 2 , 
calculated as a function of x at x' = 0, for a pair of hard point-like bosons (circles) 
and a pair of non-interacting fermions (squares) moving in ID under harmonic 
confinement. The full line gives the exact analytical result for the case of fermions. 
The error bars in the Monte Carlo data are within the size of the symbols. 
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Fig. 2. Monte Carlo results for the momentum distribution of a pair of hard 
point-like bosons (circles) and a pair of non-interacting fermions (squares) moving 
in ID under harmonic confinement. The full line gives the exact analytical result 
for the case of fermions. The error bars in the Monte Carlo data are within the size 
of the symbols. 
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Fig. 3. Log-log plot of the momentum distribution of a pair of hard point-like bosons 
moving in ID under harmonic confinement, from Monte Carlo integration (circles 
with error bars). The full line shows the asymptotic law given in Eq. (7). 



